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The quantum tunneling effects between two silver plates are studied using the time dependent 
density functional theory. Our results show that the tunneling depends mainly on the separation and 
the initial local field of the interstice between plates. The smaller separation and larger local field, 
the easier the electrons tunnels through the interstice. Our numerical calculation shows that when 
the separation is smaller than 0.6 nm the quantum tunneling dramatically reduces the enhancing 
ability of interstice between nanoparticles. 
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Metal nano-gaps offering strong surface plasmon cou- 
plings have very rich physical properties. The related 
studies have been very hot topics in the field of plamon- 
ics, e.g., single molecule surface-enhanced Raman spec- 
troscopy [JB, optical nano-antennas [3], high-harmonic 
generation [4j. The electromagnetic (EM) enhancement 
near the metal surface, which is caused by the resonant 
excitation of surface plasmon ||, is the dominating rea- 
son for the surface-enhanced Raman scattering (SERS) 
0, Q. Huge SERS with single molecule sensitivity can 
be obtained when molecules are located in the nano-gap 
between two metallic nano-structures [H 0, 0, B 0| • A 
lot of efforts have been made to seek extreme sensitive 
SERS substrates [13, HH, HI] • 

Theoretically, people have used m any methods based 
on the classical electrodynamics [l3|, 0, UH to estimate 
the SERS enhancement. These classical results indi- 
cate that the smaller the nano-gap, the higher the en- 
hancement. However, as the separation decreases to 1 
nm, the displacive current would partly become electron 
tunneling current which can reduce the EM enhance- 
ment substantially 16]. A recent experiment on the 
four-wave mixing at coupled gold nanoparticles clearly 
demonstrated that the quantum tunneling (QT) effect 
becomes significant for the distance smaller than 0.2 nm 
[]~7l ] . and a recent study of the plasmon resonance of a 
nanoparticle dimer gave quantum description of such a 
phenomenon [l^]. It is well known that the EM enhance- 
ment is the main contribution to SERS. Its enhancement 
factor is proportional to the fourth power of the local field 
enhancement, i.e. M 4 , where M=\Ei oc \/\E \ with E\ oc 
and Eq being the local enhanced electric field and the 
incident electric field, respectively. Therefore, even for 
small QT effects on M, after a fourth power, the influence 
to SERS could be huge. In this Letter we investigate the 
effects of QT on SERS with the time dependent density 
functional theory [19(. Our studies are able to quantify 
these effects and point out at exactly what conditions the 
QT has to be taken into account. 

As the "hot spot", where the SERS is strongest, is lo- 
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FIG. 1: Schematic drawing of the "hot spot" between two 
silver nano-spheres. As the "hot spot" (shaded area) is small, 
its local field is almost identical to the one computed by re- 
placing the sphere with a plate. Eq is the incident laser field. 



calized in a very small volume in the interstice between 
particles, it is convenient to investigate the QT effect 
between two closely placed plates instead of two nano- 
particles. As shown in Fig. 1, in the vicinity of the "hot 
spot" (shaded area), two plates are not much different 
from two nano-spheres. Besides we use two approxima- 
tions for our numerical calculations: (1) In the general- 
ized Mie theory, the electric current inside nano-sphere 
is set to be zero [Hj], so we can regard the silver plates 
as equipotential bodies at all time in our calculation; (2) 
The laser field is treated as a static electric field, and 
the QT effect in an oscillating field can be described by 
the results of static field in one period of laser. With 
these simplifications, when the separation d is not very 
small, the electric tunneling effect can be studied by the 
method developed by Simmons [2fj| . which regards elec- 
trons are tunneling through a voltage barrier. We find 
that Simmons' method is not proper when the distance 
d < 1 nm. For example, at d = 0.6 nm, the mean barrier 
height becomes negative at low voltage limit, indicating 
the failure of this method. 

In this work we adopt a more sophisticated method, 
the time dependent density functional theory (TDDFT) 
fl9l ] with the jellium model, where the ionic lattice is 
treated as a uniform positive charge background. In this 
method, we solve self-consistently a set of time dependent 
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Schrodinger equations, 
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^-g^ + Veff(x i t)+V ext (x i t) (1) 



where we have used the atomic units and ipk denotes a 
quantum state inside the Fermi surface of the silver plate. 
V ex t(x,t) is the external potential coming from the laser 
field and its induced field. V e ff(x,t) is the effective po- 
tential felt by an electron through Coulomb interaction 
and correlation and exchange; it depends on the elec- 
tron density. In our approach, we use Crank-Nicholson 
method [2l| to update the wave function. To quantify the 
QT effects on the SERS, we monitor time evolution of the 
potential difference SV between the two silver plates. We 
compute SV with the formula 
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where x r l and x r r , are coordinates of the left and right 
surfaces of the right (left) plate. 

Let us now turn on the laser field. The electrons inside 
each silver plate will start moving instantly to counter- 
balance the applied electric field so that the total electric 
field inside each plate is zero. At the same time, an en- 
hanced field is induced in the "hot spot" . Afterwards, 
the electrons will start to tunnel between the two silver 
plates under the following external potential 

Eqx x < X\ , 

V ext (x, t) = { E loc (x - X x ) + E X X X^x^X 2 , 
E (x - d) + Ei oc d x>X 2 . 

(3) 

It is clear from the above analysis that the initial elec- 
tron state for the Schrodinger equations in Eq.([T]) is the 
state where the electrons have moved to counter-balance 
the incident laser field. To obtain this initial state, we 
compute with the method developed by Schulte [22j the 
ground state of the metallic plate under the following 
external potential 



V°(x) = 




x < x l 

x l < x < x r 

x > x r 
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where x , x r are the left and right surfaces of the plate. 

Figure [3] shows the calculated time evolution of the 
potential difference SV between the two plates separated 
by d = 0.3 ~ 1 nm. The strength of the incident electric 
field is Eq = 2.74 x 10 5 V/m, corresponding to a laser 
with power P — 100 fj,W and focal spot ~ 1 /im. In most 
SERS experiments, even for single molecule detection, a 
much smaller P ~ 1 /iW is used P, The diameter 
of nano-particle is D = 6 nm. Note that we have calcu- 
lated for three different diameters D = 4, 5, 6 nm and 
the results are almost identical . This indicates that the 
physical process in the "hot spot" is not sensitive to ge- 
ometric features that are far away, further justifying our 
replacement of the spheres with the plates. 

We see in Fig. f2]that SV decays while oscillating with a 
frequency close to the bulk plasma frequency. The decay 



gets severe as the separation becomes smaller. This kind 
of decay can be intuitively understood by viewing the 
system as a bad capacitor that leaks current. 
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FIG. 2: Time evolutions of the potential difference 8V 
between the silver plates for different separations d = 
0.3, 0.4, 0.5, 0.6 and 1 nm. D = 6 nm; E = 2.74 x 10 5 
V/m; M = 1000. Dashed lines are for the averaged 5V. 
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FIG. 3: The decay rate rj of SV as a function of the separation 
d. D = 6 nm; E = 2.74 x 10 5 V/m; M = 1000; A = 500 nm. 

To measure the decay, or the QT suppression of the 
enhanced field, we introduce a decay rate defined by 



77 = 1- 



5V(0) 



(5) 



where T is the typical optical period of the incident laser, 
e.g., T — 1.67 fs for a laser wavelength A = 500 nm. Note 
that SV(-j) is not the value of SV but the averaged value 
of SV over one oscillation period at t = T/2. As the 
distance d decreases, the potential difference decays with 
time dramatically. At d — 0.3 nm, the local field is re- 
duced by ~ 86.6% after half optical period (A = 500 
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nm), that is, the SERS enhancement is 3.1 x 10 3 times 
((1 — ?/)~ 4 ) smaller than the one obtained from classi- 
cal theory. By contrast, at d = 1 nm, the reduction of 
the local field enhancement by the QT is only 14%, corre- 
sponding to one time decrease of SERS enhancement fac- 
tor. This means that the enhancement can be sustained 
if the separation is larger than 1 nm. The decay rates r\ 
for these different separations are computed and plotted 
in FigO where we see r\ decreases exponentially as d in- 
creases. Specifically, when the separation is smaller than 
0.6 nm, the QT can reduce the local field significantly. 
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FIG. 4: (color online) Time evolutions of the potential differ- 
ence SV between the silver plates for different enhanced local 
electric fields Ei oc . d — 0.6 nm; D = 6 nm. The inset shows 
the decay rates r\ as a function of the local electric field Ei oc . 
d = 0.6 nm; D = 6 nm. 

It is evident that both the enhancement factor M and 
the laser power P can affect the QT via the enhanced 
local field Ei oc , which is proportional to M\fP. We find 
through numerical calculations that for the range of laser 
power commonly used in experiment, the deciding factor 
is My/P, not individual values of M and P. For example, 
we find that the time evolution of 6V for P = 10 mW, 
M = 100 and P = 100 fj,W, M = 1000 is almost the same 



(not shown). This means that we need to consider only 
the enhanced local field E[ oc . Figure Q] shows the time 
evolutions of 5V for different Ei oc at d = 0.6 nm and 
D = 6 nm. We see clearly that larger local field induce 
larger QT, which in turn reduces the enhancement. As 
shown in the inset in Fig. [4] the decay rate rj decreases 
slowly when Ei oc < 2 x 10 9 V/m, and reaches a non-zero 
constant when Ei oc goes to 0. This can be explained by 
the fact that when the tunneling is small, we still have 
the linear current-voltage relation [2(J, J(t) — (35V(t). 
From this relation, we obtain 

SV(t) = 8V{Q)e~ dpt . (6) 

Therefore, when SV(0) — » 0, we have the minimal decay 
rate r\ = 1 - 5V(T/2)/5V(0) -► 1 - e ~ df}T / 2 . It should 
be noted that the minimum decay rate is determined by 
the separation. At d = 0.6 nm, the minimum reduction 
is about 26%. 

We emphasize that the reduced SERS calculated by us 
is not necessarily the overall SERS of a molecule placed 
in the nano-gap. With a molecule in the gap, the sit- 
uation can become much more complex. On the one 
hand, the oscillatory tunneling current can be coupled 
to the molecule inelastically, generating additional Ra- 
man signals 0, [24j. On the other hand, the chemical 
enhancement can also be affected by the QT [25]. Thus 
the reduced EM enhancement might be compensated or 
even over- compensated by these two factors. More stud- 
ies are needed to clarify the issue. 

In sum, we have investigated the time evolution of QT 
between two plates using the TDDFT method. We have 
found that smaller separation and larger local field result 
in stronger QT. Our numrical results show that when the 
separation is smaller than 0.6nm, the suppression of the 
EM part of enhancement in SERS is very significant for 
the common laser power used in experiment. 
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